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A MATHEMATICAL  THEORY  FOR  VARIABLE-COEFFICIENT 
LAN  CHESTER-TYPE  EQUATIONS  OF  'MODERN  WARFARE' 


by 

James  G.  Taylor  and  Gerald  G.  “Brown 


1.  Introduction. 

Deterministic  (in  the  sense  that  the  outcome  is  determined  with  certainty  by 
Initial  conditions)  Lanchester-type  equations  of  warfare  (see  [37],  [A3])  are  of  value 
for  identifying  trends  in  weapon  system  analysis  or  force  structuring  studies  because 
of  their  computational  convenience,  even  though  combat  between  two  opposing  military 
forces  is  Indeed  a complex  random  process  (see  Note  1).  The  work  of  S.  Bonder  [6], 

[8],  [11]  on  methodology  for  the  evaluation  of  military  systems  (in  particular,  mobile 
systems  such  as  tanks,  mechanized  Infantry  combat  vehicles,  etc.)  provides  a motivation 
for  interest  in  variable-coefficient  deterministic  Lanchester-type  combat  formulations 
and  having  analytic  solutions  to  these  available  (see  Note  2)  (especially  as  this 
facilitates  parametric  analysis  (see  [8])).  Bonder  [7]  has  pointed  out  that  (at  least 
for  the  case  of  mobile  weapon  systems)  the  validity  of  the  assumption  of  constant 
attrition-rate  coefficients  is  open  to  question.  [Indeed,  Bonder's  work  [6],  [8]  has 
stimulated  the  interest  of  the  authors  on  this  subject.]  Additionally,  analytic  results 
for  such  variable-coefficient  formulations  have  also  proven  useful  in  the  study  of  the 
optimal  control  of  deterministic  Lanchester-type  attrition  processes  [36],  [38],  [39]. 
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In  several  previous  papers  one  of  the  authors  has  given  analytic  solutions  for 
a few  special  cases  of  variable-coefficient  Lanchester-type  "square-law"  attrition 
equations  for  combat  between  two  homogeneous  forces.  In  [34)  (in  which  Taylor  inadvert- 
ently rediscovered  some  results  apparently  first  observed  by  B.  0.  Koopman  [35]  (see 
pp.  65-67  of  [29]))  it  was  shown  that  the  solution  to  such  equations  is  no  more  compli- 
cated than  that  for  constant  coefficients  when  the  ratio  of  attrition-rate  coefficients 
is  constant  (see  also  [11]).  In  [37]  more  general  results  were  given  via  successive 
approximations  and  applied  to  power  attrition-rate  coefficients  (both  for  a case  which 
may  be  thought  to  apply  for  two  weapon  systems  with  the  same  maximum  effective  range 
and  for  one  with  coefficients  linear  in  time  but  reflecting,  for  example,  different 
effective  ranges).  The  infinite  series  solutions  given  in  [37]  only  apply,  moreover, 
to  th?  restrictive  case  in  which  the  lethality  of  at  least  one  side's  fire  is  initially 
zero.  One  interpretation  of  this  is  that  the  battle  must  begin  at  the  minimum  of  the 
maximum  effective  ranges  of  the  two  systems.  The  new  theoretical  results  of  the  paper 
at  hand  allow  these  previous  results  to  be  extended  to  more  general  cases  of  interest. 

Consideration  is  given  to  solution  in  terms  of  tabulated  functions.  Since 
tabulations  of  such  functions  tor  the  cases  at  hand  do  not  exist,  we  introduce  Lan- 
chester  functions  that  should  prove  useful  in  this  respect  and  give  their  mathematical 
properties  that  facilitate  solution.  Additionally,  we  determine  what  assumptions  about 
the  attrition-rate  coefficients  lead  to  simplification  of  analytic  results. 

The  organization  of  this  paper  is  as  follows.  We  first  consider  the  general 
model  for  which  analytic  results  are  developed  and  then  discuss  the  concept  of  solution 
in  terms  of  tabulated  functions.  Hie  main  results  of  this  paper  are  the  development  of 
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a mathematical  theory  of  Lanches ter- type  equations  for  a "square-law"  attrition  process 
for  combat  between  two  homogeneous  forces.  These  general  results  are  then  applied  to 
some  special  types  of  attrition-rate  coefficients  of  interest,  and  new  mathematical 
functions,  Lanchester  functions,  which  could  be  tabulated  in  the  future  are  proposed. 
Some  numerical  examples  are  given.  Then,  we  show  how  the  classical  Lanchester  square 
law  generalizes  to  such  variable-coefficient  formulations.  Finally,  we  discuss  exten- 
sions, the  significance,  and  applications  of  our  results. 


2.  Lanchester-Type  Equations  of  Modern  Warfare. 


In  1914  F.  W.  Lanchester,  an  English  aeronautical  engineer  who  lived  from  1868 
to  1946  (see  [27]),  in  order  to  provide  insight  into  the  dynamics  of  combat  under 
"modem  conditions"  and  justify  the  principle  of  concentration  (see  Note  3),  hypothe- 
sized that  combat  between  two  opposing  forces  could  be  modelled  by  [25]  (see  Note  4) 


dx/dt  n -a(t)y 
dy/dt  * -b(t)x 


with  x(t-O)  - xn, 

° (1) 
with  y(t-O)  - yQ, 


where  x(t)  and  y(t)  denote  the  numbers  of  X and  Y combatants  and  a(t)  and 
b(t)  are  (today)  called  Lanchester  attrition-rate  coefficients.  These  coefficients 
represent  the  lethality  of  each  side's  fire. 
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We  shall  refer  to  the  equations  (1)  as  Lanchester's  equations  of  modem  warfare. 


IWo  sets  of  circumstances  under  which  these  equations  have  been  hypothesized  to  apply 
are  as  follows: 

(a)  both  sides  use  aimed  fire  and  target  acquisition  times  are  constant,  Independent 
of  the  force  levels  (a  special  case  is  when  they  are  negligible)  [43], 

(b)  both  sides  use  area  fire  and  a constant  density  defense  [12]. 

A more  complete  discussion  of  these  hypotheses  (for  constant  lethality  of  fires)  is  to 
be  found  in  the  above  referenced  original  papers.  Other  factors  may  be  included  in 
the  equations  and  other  diffeiential  equation  models  of  combat  attrition  may  be  referred 
to  as  Lanches ter- type  equations,  but  we  will  not  consider  these  here  (see  [16],  [37]). 

In  Lanchester's  original  work  [25],  the  attrition-rate  coefficients  were  assumed 
to  be  constant.  When  either  this  is  true  or  their  ratio  is  constant,  i.e.  a(t)/b(t)  ■ 
kg/k^,  the  classical  Lanchester  square  law  results 

kb(x2-x2(t))  = ka(y2-y2(t)).  (2) 

which  has  the  important  implication  that  a side  can  significantly  reduce  its  casualties 
by  initially  committing  more  forces  to  battle.  As  w°  show  below  by  consideration  of 
the  force-ratio  equation,  a generalization  of  (2)  holds  even  for  the  general  case  of 
variable  coefficients  in  which  a(t)/b(t)  is  not  constant.  Thus,  it  seems  appropriate 
to  refer  to  (1)  as  the  equations  for  a "square-law”  attrition  process  (see  also  [36], 
[39]). 

TVo  significant  developments  in  the  Lanchester  theory  of  combat  during  th ».  1960's 
were  (a)  the  development  of  methodology  for  the  prediction  of  Lanchester  attrition-rate 
coefficients  from  weapon  system  performance  data  by  S.  Bonder  [7],  [9]  and  others  [3], 
[11],  and  (b)  G.  Clark's  development  of  methodology  for  the  (maximum  likelihood)  estima- 
tion of  such  coefficients  from  Monte  Carlo  simulation  data  [13].  Both  these  developments 
(and  others  [31],  [32])  have  facilitated  the  application  for  defense  planning  studies 
of  models  such  as  (1)  and  its  generalizations  to  combat  between  heterogeneous  forces 
(see  [11]). 
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Very  recently,  F.  Grubbs  and  J.  Shuford  [19]  have  applied  concepts  from  reliability 
theory  to  develop  a new  probabilistic  formulation  for  Lanchester  combat  theory.  They 
consider  the  random  combat  process  to  be  modelled  by  the  expected  fraction  of  survivors, 
i.e.  x/xq  and  y/y^,  being  functionally  related  to  various  forms  of  probability 
distributions  for  the  random  time  to  kill.  Although  they  consider  mainly  cases  in  which 
friendly  losses  are  not  related  to  the  enemy  force  level,  they  do  set  up  the  following 
model  (which  is  not  solved) 

d[(x0-x)/x0J/dt  - f (t)y/yQ  and  d[(y0-y)/yQ j/dc  - g(t)x/xQ,  (3) 

where  f(t)  and  g(t)  are  the  time-to-kill  probability  density  functions.  They  suggest 
that  the  Weibull  distribution  is  a convenient  point  of  departure  for  the  latter.  Our 
general  results  here  may  be  applied  to  the  formulation  (3)  suggested  by  these  authors 
[19].  It  should  be  pointed  out  that  a significant  advantage  of  the  fomulation  (3)  over 
the  usual  continuous  parameter  Markov  chain  formulation  is  one  of  computational  con- 
venience (see  [13]). 

The  Lanchester  attrition-rate  coefficients  in  (1)  depend  on  a number  of  variables 
such  as  firing  doctrine,  firing  rate,  rate  of  target  acquisition,  force  separation,  ’ 
tactical  posture  of  targets,  etc.  (see  [6]  or  pp.  18-26  and  pp.  81-114  of  r 11  ] for 
attrition-rate  coefficient  prediction  methodology).  Bonder  [6]  (see  also  [11])  has 
considered  a number  of  forms  for  attrition-rate  coefficients  baqed  on  examination  of 
data  for  some  representative  weapon  systems.  Motivated  by  this  work,  we  will  consider 
the  following  coefficients  in  the  paper  at  hand 

(I)  a(t)  * ka(t+C)m  and  b(t)  ■ k^(t+C)n  with  C 0,  (4) 

(II)  a(t)  - ka(t+C)  and  b(t)  - k^t+C+A)  with  A,C  * 0.  (5) 

Some  situations  which  may  be  modelled  with  these  coefficients  are  discussed  below. 
[Additionally,  for  m and  n negative  (for  m,n  < -1  we  require  C > 0)  equations 
(1)  with  (4)  may  be  considered  to  model,  for  example,  an  infantry  "fire  fight’’  in  which 
the  combatants  "take  cover"  so  that  the  lethality  of  fires  decreases  with  time.] 


5 


The  results  given  In  [37]  were  for  the  case  In  which  C ■ 0 (sec  equations 
(4)  and  (5)  above).  For  C > 0 the  methods  of  [37]  (both  successive  approximations 
and  Infinite  series)  have  not  been  mathematically  tractable.  However,  the  theoretical 
results  of  the  paper  at  hand  considerably  simplify  the  situation. 

3.  Tabulated  Functions  Arising  in  Military  Operations  Research. 

It  seems  appropriate  to  discuss  two  points:  (a)  the  Importance  of  solution 

in  terms  of  "tabulated  functions,"  and  (b)  why  (1)  with  cither  (4)  or  (5)  does  not 
lead  (in  general)  to  previously  tabulated  functions.  This  should  help  provide  perspec- 
tive on  the  contributions  of  the  paper  at  hand. 

Let  us  consider  the  case  in  which  a(t)  and  b(t)  ate  constants,  e.g. 
a(t)  « a * constant.  Then,  we  say  that  we  have  "solved"  for  the  X force  level  when 
we  obtain  x(t)  ■ cosh»ab  t - y^/a/b  sinh»^b  t.  This  result  is  useful  because  the 

solution  is  expressed  in  a form  (i.e.  in  terms  of  two  functions  of  a single  argument 
and  tabulations  of  these  functions  are  available  (see,  for  example,  [1]))  that  allows 
u8  to  conveniently  generate  numerical  results.  Thus,  we  can  generate  a numerical 
value  for  the  X force  level  for  given  values  of  the  parameters  x^,  y^,  a,  b,  and 
t by  computing  the  value  of  the  argument,  i.e.  *^ab  t,  extracting  (possibly  using 
interpolation)  the  corresponding  values  of  the  hyperbolic  functions  from  tables,  and 
some  simple  algebraic  manipulation.  [Although  a representation  of  such  a transcendental 
function  as  an  infinite  series  may  be  useful  for  generating  such  tables,  it  is  clear 
that  the  availability  of  such  tables  is  highly  desirable.]  One  goal  of  our  research 
has  been  to  reduce  the  solution  of  (1)  (with  coefficients  like  (4)  or  (5))  to  as  close 
as  possible  to  the  above  situation. 

in  trying  to  develop  a solution  for  (1)  with  (4)  or  (5),  the  authors  have,  in 
general,  not  been  able  to  express  results  in  terms  of  previously  tabulated  functions 
(see  below).  This  is  not  surprising  when  one  recalls  that  a significant  number  of  the 
so  called  "special  functions"  (sec  [1])  arise  in  the  solution  of  equations  of 
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mathematical  physics  In  various  coordinate  systems  (see  Chapter  5 of  [28]).  Moreover, 
the  X force-level  equation  (6)  with  coefficients  (5)  appears  to  be  one  not  previously 
encountered  (see  Note  S). 

Thus,  one  contribution  of  this  paper  Is  to  point  out  that  essentially  new  special 
functions  are  Involved  In  solving  the  above  variable-coefficient  Lanches ter- type 
equations  and  that  tabulations  of  the  Lanchester  functions  introduced  below  would 
facilitate  this.  It  should  be  pointed  out  that  a number  of  "special  functions"  and 
various  tabulations  of  these  have  previously  arisen  in  target  coverage  problems  of 
military  operations  research  (see,  for  example,  [17]). 


4.  A Mathematical  Theory  of  Lanches ter-Type  Equations  of  Modem  Warfare. 

The  Lanchester-type  equations  (1)  yield  the  X force-level  equation 


with  initial  conditions 

i(t-O)  - x0.  and  fjTo  ■ ‘V 

The  solution  of  (2)  is  given  by 


(6) 


x (t ) - CjX^t)  + C2x2(t),  (7) 

where  (x^(t) ,x2(t) } denotes  a fundamental  system  of  solutions  (see  p.  119  of  [23]). 
These  functions  can  be  chosen  tc  satisfy 

dx  jv, 

“ ka(t)y2,  and  ~ - - ka(t)ylt  (8) 

where  (yj(t),y2(t))  denotes  a fundamental  system  of  solutions  to  the  Y force-level 
equation.  The  complementary  Y-functions  y^(t)  and  y2(t)  consequently  satisfy 
equations  similar  to  (8). 
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In  (8)  above,  the  constant  k la  completely  arbitrary  and  depends  on  how  the 
fundamental  X-  and  Y-f unctions  are  defined.  It  is  convenient  to  define  these  so 
that  they  reduce  to  well-known  elementary  (transcendental)  functions  In  special  cases. 
One  such  special  case  Is  when 


a(t)  - kah(t),  and  b(t)  - l^hCt).  (9) 

In  this  case  it  is  natural  to  take,  for  example, 

XjCt)  ■ cosh  fl(t),  and  y2(t)  ■ 8inh  0(t),  (10) 

where  0(t)  ■ v^a^b  \ h(s)ds  + 0(t«O).  It  follows  that  k ■ A^/kfl. 

The  complementary  X-  and  Y-functions  possess  the  rather  remarkable  property 
that  (Independently  of  the  determination  of  k) 


xl(t)yi(t)  - x2(t)y2(t)  - 1 Vt. 


(11) 


Hits  propercy  (11)  is  essential  for  determining  the  constants  and  C2  of  (3)  in 

the  general  case  in  which  x^(t-O)  r 0 and  x2(t**0)  i 0.  The  p-.^pertles  of  the 
general  Larchester  functions  are  summarized  in  Table  I.  In  general,  the  function  x^(t) 
and  y^(t)  are  sort  of  like  the  hyperbolic  cosine  (with  the  appropriate  argument), 
while  x2(t)  and  y2(t)  are  sort  of  like  the  hyperbolic  sine  (see  (13)  below). 

Let  us  now  sketch  the  proof  of  (11).  By  Abel's  identity  (see  p.  75  of  (23]), 
the  Wronskian  of  Xj,  denoted  as  W(x^,x2),  satisfies 


W(xltx2) 


X1  X2 


a ex 


P<-  | l aTtT  ’ 


whence  for  all  t 


a ■ k{x1(t)y1(t)  - x2(t)y2(t)},  (12) 

where  a Is  a constant  to  be  subsequently  determined.  It  will  be  convenient  to  choose 
the  above  general  Lanchcster  functions  so  that  at  t * t^  we  have 

xl(t"t0)  " y^t-tQ)  = 1,  and  x2(t-t0>  - y2(t-t0)  - 0.  (13) 

It  follows  that  a ■ k,  and  the  proof  of  (11)  is  complete. 
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TABLE  I 


Properties  of  the  Ceneral  Lenchester  Functions 


b(t)x2 


4. 


b(t)x1 


The  general  Lanchestcr  functions  , x^i  and  y,  my  be  determined  by 

either  successive  approximations  or  the  method  of  Frobenious  (see  [37]).  The  time  tg 
is  chosen  so  that  the  general  Lanchester  fi  actions  take  a convenient  form.  Examples 
of  this  are  given  below. 

The  constants  and  in  (7)  are  determined  by  the  Initial  conditions  for 

(6),  and  thus  via  (8)  and  (11)  it  follows  that 

x(t)  - x0{y1(t-0)x1(t)-y2(t-0)x2(t)} 

- y0*ir7k^(xi(t-0)x2(t)-x2(t=0)xi(t) }.  (14) 

The  expression  (14)  is  the  generalization  of  the  well-known  form  of  the  X force  level, 
x(t),  for  the  classical  Lanchester  (constant  coefficient)  equations  of  "modern  warfare" 
to  the  most  general  case  of  variable  Lanchester  attrition-rate  coefficients.  The 
reader  should  note  the  key  role  that  (11)  plays  in  verifying  that  (14)  does  indeed 
satisfy  the  initial  conditions  to  (6).  Ve  also  see  from  (14)  thal  once  results  are 
obtained  for  one  member  of  a family  of  battles,  force  levels  may  be  readily  deduced 
for  other  members  of  this  family. 

Recalling  the  well-known  constant  coefficient  results,  the  reader  knows  that 
for  tg  + 0 (14)  is  at  least  sometimes  capable  of  further  simplification.  However, 

it  is  indeed  remarkable  that  this  is  only  true  when  (9)  holds. 

THEOREM  1:  For  tg  i 0,  any  further  simplification  ir  (14)  is 

possible  if  and  only  if  (9)  holds  (constant  ratio  of 
attrition-rate  coefficients). 


PROOF:  From  (8)  we  have,  for  example,  that 

dx 


_L  . V2  siu  yJ. 

dy-  b(t)  x. 


Thus  a relationship  that  is  independent  of  t exist:;  between  Xj  and  y 2 if  and  only 
if  a(t)/b(t)  ■ constant.  Moreover,  when  such  a relationship  does  exist,  we  have  via 


(9)  that 
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(15) 


/ (x^t)]2  - fy2(t))2  - 1. 

Besides  being  a classical  Lanchester  "square  law,"  equation  (IS)  is  a necessary  and 
sufficient  conditio.',  for  example,  for  x^(t)  to  posseas  a so-called  algebraic  addi- 
tion theorem  (see  Note  6).  Unless  the  general  Lanchesto-  functions  possess  such 
algebraic  addition  theorems,  for  tg  j*  0 there  is  no  further  simplification  to  (14). 

Q.E.D. 

Remark  1:  When  (9)  does  hold,  we  have  that  x^(t)  ■ y^(t)  and  x2(t)  ■ y2(t),  where 

x^(t)  and  y2(t)  are  given  by  (10),  so  that  the  complementary  X-  and  Y-functions 
are  (only  in  this  case)  identical.  Moreover,  (11)  (or,  equivalently,  (15)  is  then  a 
wel.  -known  property  of  the  hyperbolic  functions.  Furthermore,  in  this  case  (14)  becomes 
the  well-known  result  [11],  [29],  [34],  [35] 

t t 

x(t)  - Xq  cosh(^Tk^  J h(s)ds)  - yQ •^ca/kb  sinh  J h(s)ds),  (16) 

0 0 

thanks  to  the  well-known  (algebraic)  addition  theorems  for  the  hyperbolic  functions, 
for  example,  like  cosh  u cosh  v - sinh  u sinh  v ■ cosh  (u-v).  It  is  also  enlightening 
to  observe  that  (16)  may  be  written  as 

x(t)  - Xq  cosh(«4(t)b(t)  t)  - y0/a(t)/b'(t)  sinh(>4(t)b(t)  t),  (17) 

__  . l rt 

where  /a(t)b(t)  denotes  an  average,  i.e.  /a (t)b(t)  * — /a(s)b(s)  ds.  Thus,  we 

C J0 

see  that  in  the  case  of  a constant  ratio  of  attrition-rate  coefficients,  the  X f rce 
level,  x(t),  is  given  by  an  expression  formally  equivalent  to  that  for  the  constant 
coefficient  case  with  averages  being  used. 

Remark  2:  In  cases  in  which  a (square-law)  relationship  like  (15)  does  not  hold,  there 

is  no  such  convenient  reduction  of  (14)  with  tQ  i 0 to  a simpler  form  like  (16)  via 
an  algebraic  addition  theorem.  Thus,  in  general  the  X force  level,  x(t),  does  not 
take  a simple  form  except  when  tg  ■ 0 or  a(t)/b(t)  ■ constant. 
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Remark  3i  When  tg  • 0,  then  (14)  again  takes  a particularly  simple  lorn 

x(t)  - *0*1<t)  “ y^kjk^  x2(t).  (17) 

Our  previous  results  [37]  were  all  of  the  form  (17). 

In  a previous  paper  [37]  (see  also  [34],  [35])  we  showed  that  when  the  ratio 
of  attrition-rate  coefficients  is  constant,  we  can  transform  the  • v force-level  equa- 
tion into  one  with  constant  coefficients  by  a transformation  of  the  independent  variable 
t.  As  we  have  seen,  this  case  leads  to  particularly  convenient  results.  In  this  respect, 
a useful  theorem  is 


THEOREM  2:  A necessary  and  sufficient  condition  to  be  able  to  transform 

the  X force-level  equation  (6)  by  a transformation  of  the 
Independent  variable  t into  a linear  second  order  ordinary 
differential  equation  with  constant  coefficients  is  that 


1 db 


la(t)  dt  bit)  dt 


} “ CONSTANT. 


PROOF:  The  theorem  follows  immediately  from  a result  given  on  pp.  73-74  of  [4].  Q.E.D. 

Moreover,  when  the  X force-level  equation  is  so  transformable,  the  desired  sutstitu- 
tion  is  given  by  u » K J /a(s)b(s)  ds,  where  j ...  ds  denotes  an  indefinite 
integral  and  K is  an  arbitrary  constant  conveniently  chosen. 

For  example,  for  the  power  attrition  rate  coefficients  (4),  Theorem  2 tells  us 
that  we  can  transform  (6)  into  an  equation  with  constant  coefficients  only  when  (I) 
m » n,  or  (II)  m + n + 2 - 0. 

5.  Power  Attrition-Rate  Coefficients  With  Both  Systems*  Effectiveness  Zero  at  the 
Same  Time. 

In  this  case  the  attrition-rate  coefficients  arc  given  by  (4).  There  arc  two 
cases  to  be  considered,  depending  on  whether  (6)  with  (4)  can  be  transformed  into  a 
constant  coefficient  equation  (see  Theorem  2 above): 
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(I)  m + n + 2 i 0, 
and  (II)  b + n + 2 ■ 0, 


Case  I.  m + n + 2 0: 

In  order  for  solution  methods  (either  successive  approximations  or  the  method 
of  Frobenlus)  to  be  applicable,  we  must  further  Impose  the  following  restrictions: 

(a)  for  C ■ 0,  we  must  have  m > -1  and  n > -1,  while  (b)  for  C > 0,  we  must 
have  m + n + 2i*0.  It  should  be  noted  that  only  the  former  case  was  considered  In 
our  previous  paper  [37]  and  that  the  theory  presented  in  the  last  section  Is  essential 
for  extending  these  results  to  the  latter  case. 

From  results  given  In  [37],  it  follows  that  a fundamental  system  of  solutions 
to  (6)  with  attrition-rate  coefficients  (4)  is  given  by 

xx(t)  - r(q)(/Tk^/2s)l-'(t+cJm+1)/2  I_p(4Tk^(t+C)8/s),  (18) 

and 

x2(t)  - r(p)(/lTk^/2s)q(t+cJ0+1)/2  IF(/n£(t+C)8/s),  (19) 

where  I denotes  the  modified  Bessel  function  of  the  first  kind  of  order  p, 

P 

p ■ (m+l)/(mfn+2) , p + q ■ 1,  and  2s  ■ m + n + 2,  with  similar  results  holoing 
for  ly1(t),y2(t)}.  The  above  Lanchester  functions,  of  course,  satisfy  all  the 
properties  given  in  Table  I The  X force  level,  x(t),  then  is  given  by(14). 

To  emphasize  the  dependence  on  the  parameter  C it  is  sometimes  convenient  to 

write,  for  example,  x^(t)  ■ x^(t;C).  It  should  bf  noted  then  chat  Xj(t*0;O0)  ■ 
y^(t"0;O0)  i and  x2(t*0;O0)  ■ y2(t-0;C=0)  ■ 0 so  that  (14)  with  these  particular 
fundamental  systems  reduces  to  (17)  when  C ■ 0.  It  should  be  noted  that  no  such 
simplification  of  (14)  occurs  for  C > 0 unless  m ■ n.  In  other  words,  the  param- 
eter tj  with  the  property  (13)  is  given  by  tg  - -C. 

The  above  representations  (18)  and  (19)  of  fundamental  systems  of  solutions  to 

the  X and  Y force-level  equations  are  not  particularly  convenient  because,  for 
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example,  for  m,n  > -1  we  have  0 < p,q  and  p + q ■ 1,  and  tabulations  only  exist 
of  the  modified  Bessel  function  of  the  first  kind  of  fractional  order  v for  a restric- 
tive set  of  values  (l.e.  v ■ ±1/4, ±1/3, 12/3, 13/4  (see  (37))).  This  Bessel  function 
reduces  to  other  tabulated  forms  (i.e.  hyperbolic  functions)  for  v - ±1/2.  Therefore, 
it  is  more  useful  to  express  x^(t),  x^Ct),  y^(t),  and  y ( t ) in  the  form  of  infinite 
series  and  to  consider  the  resulting  transcendental  functions  as  entities  in  their 
own  right.  Thus,  it  is  convenient  to  define  the  folic iwng  Lanchester  functions 


tn,n 


m,n 


m,n 


m,n 


(t) 

■ r(q) 

•ft 

pvb' 

tk(mfn+2) 

(20) 

o 

-O  H 
JC 

m+iY+2j 

k ! r (k+q)  * 

(t) 

- r(p) 

CO 

Y 

tk(nHm+2)+m+l 

(21) 

L 

k*0 

m+n+2j 

k ! T (k+l+p) 

(t) 

- r(P) 

y 

rSi 

^ tk(m+n+2) 

(22) 

L 

k-0 

|m4-n+2j 

k!I’(k+p)  ' 

(t) 

- T(q) 

00 

Y 

\A  k ) 

a h 

2k+1 

tk(nrHi+2)+n+l 

(23) 

l 

k-0 

m+n+2j 

k ! T (k+l+q) 

The  Lanchester  functions  u , v , U , 
m,n  ra,n  m,n 


and  V have  tl.c*  properties 
m , n 


shown  in  Table  II.  For  reasons  to  be  explained  below  we  will  refer  to  u^  ^(t)  and 

U (t)  as  complementary  Lanchester  functions  and  similarly  for  v (t)  and 
m,n  m,n 


V (t).  The  solution  (14)  to  (6)  may  then  be  written  in  terms  of  these  Lanchester 

tn,n 


functions  as 


x(t)  - xn(U  (C)u  (t+C)-V  (C)v  ( t+C) ) 
U m,p  m,n  m , n m , n 


- y./k  / k.  (u  (C)v  (t+C)-u  (t+C)v  (C)),  (24) 

at)  m,n  m,n  m,n  m,n 


and  similarly  for  >(t).  From  (24)  we  see  that  the  methods  of  (37)  would  become 
hopelessly  bogged  down  in  details  for  C > 0. 

For  computational  reasons  it  is  convenient  to  introduce  the  auxiliary  Lanchester 
functions  (also  referred  to  as  the  Lanchcster-Cl i f ford-Schldf ] 1 functions,  or  LCS 
functions  (see  Note  7)),  which  are  defined  for  v / 0,-1, -2,-3, . . . 
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Table  II.  Properties  of  the  Lanchester  Functions 

u , v , U , and  V . 

«,n  n,n  B,n  n,n 


du 


»•  “St2  ■ \\  ‘VX.n(t> 


dv 


2-  -St2  ■ ^ 


dU 


“Sf  ■ /iA  VX.A 


dV 


*•  "Sf  ■ ''‘‘A  lV°lu.,n(t) 


5. 


6. 


7. 


8. 


9. 


u 

m,n 


u 

m,n 


v 

m,n 


u 

m,m 


v 

m,m 


(t)U  n(t)  - v (t)V  ft) 
m,n  m,n  n(n 


(t*0)  - U (t*0)  - 1 
cnf  n 

(t*0)  - v (t**o)  ■ o 

ra,n 

(t)  * Um,o(t)  " "08h 

(t)  - V (t)  - sinh  A k. 

m,m  a d 


- 1 Vt 


t”^1/  (nH-1)  ) 
t^/CnH-l)) 
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(25) 


Fu(x)  “ I (x/2)2k/(k!  * (j+v)}, 


k-1 


k-0 


and 


J-o 


C (x)  ■ l (x/2)2k+1/(k!  ir  (J+v)}, 


(26) 


k-0 


N 


J-0 


where  we  have  adopted  the  convention  that  it  f - 1 for  N < M.  The  LCS  functions 

J-M  J 

possess  the  properties  shown  in  Table  III. 

The  Lanchester  functions  u , v , U . and  V may  be  expressed  in 
■ ■ m,n  m,n  m,n  tn,n 

terms  of  the  LCS  functions  as  follows 


u „(t>  " Fn(*(t>). 

m,n  q 


v ( t ) - t(m'n)/2  C (*(t)), 

m,n  p 


(27) 

(28) 


Um  (t)  - F (*(t)), 
m,n  p 


V_  ft)  - L(n~m)/2  G (ij»(t)) , 


m,n 


(29) 

(30) 


where  ^(t)  - f\\ k t^Jlk-t")  t/((m+n+2)/2).  Thus,  u (t)  and  U (t)  have  been 
a b tr/.i  m,n 

called  complementary  Lanchester  functions  because  of  the  above  relationships,  the 
fact  that  p + q » 1,  and  property  3.  of  Table  III.  Additionally,  the  introduction 
of  the  LCS  functions  sheds  light  on  the  parametric  dependence  of  solutions:  there 

are  two  exponent  parameters,  p and  (m-n)/2,  and  an  "intensity"  parameter,  0(0  ■ 
/a (t)b(t)  (t+C)/((m+n+2)/2),  as  well  as  the  relative  effectiveness  parameter 
■Vv  We  finally  note  that  (24)  may  now  be  written  as 


x(t)  - xn( F (0 (0) )F  (0 ( t) )-( 1+t /C) (m-n) 7 2G  (9 (0) )G  (0 ( t) ) } 

0 p q q P 

- y0/n(t*0)/b(t=0){(l+t/C) (m"n)/2Fq(e(0))Gp(e(t))-Cp(e(0))Fq(0(t))),  (31) 

where  0(t)  * /i"(t)b(t)  (t+C)  / ( (nr+n+2) /2)  . It  should  be  observed  that  average 
"intensity"  /a\t")b"(t)  is  related  to  0(t)  * 0 ( t ; C)  by 

/a ( L ) b ( t ) t = 0(t;C)  - 6(t=0;C). 


16 


Table  III.  Properties  of  the  LCS  Functions  F^Cx)  and  C^(x). 


df  " Gv<*> 


dGv  G (x) 

-p— - ■ F (x)  * (2v-l)— 

dx  v x 


F (x)F  (x)  - G (x)G  (x)  - 1 Vt 
P <1  P <1 

where  p + q - 1 and  p,q  r 0,-1, -2, -3,. 


Fv(x-0)  - 1 


Gv(x-0)  - 0 


V*>  i 

x 2v 


(x-0)  - 0 


dG 

~ (x-0)  - -~- 
dx  2v 


Fi/2<x)  • cosh  x 


^1/2 Cx)  - slnh  x 
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Without  an  algebraic  addition  theorem,  however,  this  does  not  lead  to  simplification 
of  (31).  We  do  find,  though,  that  for  C » 0 with  m,n  > -1  we  have  /a(t)b(t)  t ■ 
0(t;OO),  and  (31)  may  be  written  like  (17)  as 

x(t)  - xQ  Fq(/a(t)b(t)  t)-yQ»/aTo7bTtT  Gp(^(t)b(0  t) . (32) 

It  may  be  shown  that  lim  /a(t)/b(t)  *G  (/!a(t)b(t)  t)  ■ 0. 

t+0  p 

Considering  (31)  and  the  above,  there  is  an  interesting  way  (with  some  similari- 
ties with  the  constant  coefficient  situation)  of  thinking  about  combat  between  two 
homogeneous  forces  described  by  (1)  and  (4):  average  combat  "intensity"  as  well  as 

how  this  changes  over  time  determines  the  course  of  combat  besides  the  initial  force 
levels,  Xq  and  y^,  and  /a(t*0)/b(t"0) . 

Case  II.  m + n + 2 ■ 0 and  C > 0 : 

It  should  first  be  noted  that  the  solution  obtained  for  Case  1 becomes  indeter- 
minate (consider,  for  example,  what  happens  to  u (t)  as  defined  by  (20)).  To 

m,  n 

solve  (6)  with  (4)  and  m + n + 2 « 0,  we  make  the  substitution  u * in  t to  trans- 
form (6’'  into  an  equation  with  constant  coefficients.  Solving  this  equation,  v*  find 
that  one  fundamental  system  (with  the  properties  given  in  Table  I)  for  (14)  is  given  by 


and  similarly  for  ( y ^ ( t ) .yjft) } , where  0 - (m+1 ) / 2) ^ , a+  = 0 + (ra+l)/2  e 

-g  , and  a « -0  + (m+l)/2  ■ -8+.  In  this  case  (14)  becomes 
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which  may  also  be  written  In  the  more  convenient  form 


Although  the  above  solutions  appear  complex,  they  are  readily  evaluated  with  the  help 
of  a ("hand-sized")  portable  electronic  calculator  such  as  Is  commercially  available 
today.  (In  fact,  such  a calculator  can  even  be  "programmed"  to  facilitate  parametric 
analyses. ) 

6 • Weapon  Systems  With  Different  Effective  Ranges:  Linear  Attrition-Rate 

Coefficients. 

Another  situation  of  interest  in  which  to  apply  our  general  theory  is  that  of 

combat  between  two  homogeneous  forces  which  use  weapons  with  different  effective 

ranges.  Let  us  consider  the  example  previously  examined  by  us  [37] (sec  also  [11]) 

of  a constant  speed  attack  of  a mobile  force  against  a static  defense.  The  weapon 

systems  of  the  two  sides  have  different  effective  ranges,  and  the  lethality  of  each 

side's  fire  depends  linearly  upon  range.  We  assume  that  the  opening  range  of  battle, 

Rn,  is  £ minimum  (R  , R„) , where  R denotes  the  maximum  effective  range  of  the 
O a p a 

Y system.  The  Lanchester  attrition-rate  coefficients  for  such  a battle  may  be  written 

as  (5).  The  parameter  C is  related  to  the  opening  range  of  battle  in  comparison 

with  minimum  (R  ,R„),  whereas  A reflects  the  difference  in  maximum  effective  ranges, 
o p 

The  range  dependencies  of  these  Lanchester  attrition-rate  coefficients  and  the  opening 
range  of  battle  is  shown  in  Figure  1. 

From  results  given  in  [37],  it  follows  that  a fundamental  system  of  solutions 
to  (6)  (and  also  one  for  the  Y force-level  equation)  with  the  attrition-rate  coeffi- 
cients (5)  (which  has  the  properties  shown  in  "’tble  I)  is  given  by 

Xj(t)  - f(t+C),  x2(t)  = g(t+C) , yj(t)  = F(t+C) , and  y2(t)  - G(t+C),  (35) 
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20 


Figure  1.  Linear  attri;ion-rate  coefficients  for  weapon  systems  with  different  effective  ranges 
[Notes:  1.  The  maximum  effective  ranges  of  the  two  weapon  systems  are  denoted  as 
R and  R0.  2.  The  opening  range  of  battle  is  denoted  as  R_  and  (as  shown) 

0 Ip  U 

Rn  < minimum  (R  ,R0).] 
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(36) 

(3>) 

(38) 

(39) 


and  the  coefficients 
are  given  by  - 1 


and  C*6  are  given  in  [37]. 
on 

and  for  n > 0 


The  coefficients 


n 


{ l> 

dJ^  - ' l4n(4n-2)/(4n-k)][D^_1/(4n-k-2))+(Dj'J/(4n-k-l))]. 

I (4/3)[(4n-2)/(3n-l)]D”“J. 

n-l 

and  E®  - 1,  E*  - 2,  and  for  n > 0 


for 

k - 

o. 

for 

1 < 

k £ n 

for 

k - 

n; 

( 

E * - < [4n(4n+2)/(4n-k+2)]((E*  1/(4n-k))+(Ej‘"}/(4n-k+l))  ] , 
n 1 n-l  n-i 

' (^/3)[(4n+2)/(3n+l)]E"  ., 

n-i 


for  k-0, 
for  1 < k i n, 
/or  k - n + 1. 


It  should  be  noted  that  f(t-O)  - F(t-O)  - 1 and  g(t-O)  - C(t-O)  - 0. 

To  emphasize  the  dependence  on  the  parameter  A,  it  is  sometimes  convenient 

to  write,  for  example,  f(t)  - f(t;A).  It  should  be  noted  that  f(t;A-0)  - F(r;A-0)  - 

cosh  (A  k.  t2/2)  and  g(t;A-0)  - G(t;A-0)  - sinh  (/k  k t2/2).  For  A - 0,  the 
a d a d 

fundamental  property  f (t)F(t)  - g(t)G(t)  - 1 of  these  Lanchoster  functions  becomes 
a well-known  property  of  the  hyperbolic  functions. 

As  was  the  case  for  the  power  attrition-rate  coefficients  of  the  previous 
section,  for  computational  reasons  it  is  convenient  to  introduce  the  following 
auxiliary  Lanchoster  functions  (depending  on  two  parameters) 
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h(X,y) 

m 

" l 

n-0 

X2n 

(2n)l  k 

1 • 
-0  n 

(AO) 

v(X,y) 
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a r — ? | 
O 

X21*1 
(2n+l) ! 
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k-0  n 

(41) 

H(X,y) 

CO 
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n-0 

X2n 
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(42) 

m 

. 2n+l 

n+1  k k 

W(A,y) 

- 1 

n-0 

A 

(2n+l) ! 

l E y . 
k-0  n 

(43) 

These  auxiliary  Lan '.Hester  functions  possess  the  following  properties: 

1.  h(X,y)H(\,y)  - w(X,y)W(X,y)  - l VX,y 

2.  h(X,0)  - H(X,0)  - cosh  X, 

2.  v(X,0)  - W(X,0)  - sinh  X. 

Using  the  offset  linear  auxiliary  Lanchester  functions,  the  X force  level 
may  be  written  as 

X(t)  - xo{H(e(o),«(O))h(e(t),6(t))-w(e(o),6(O))w(0(t),5(t))} 

- yo^a/kb{h(0(o),«(o>>w(8(t>,«(t))-w(0(O),6(°))h(0(t),6(t))},  (44) 

where  6(t)  ■ A/(t+C)  and  0(t)  * (t+C)2/2.  The  expression  6(t)  reflects 

how  much  the  above  Lanchester  functions  deviate  from  the  hyperbolic  functions. 


Some  Numerical  Exc 


In  this  section  we  examine  three  numerical  examples  which  illustrate  possible 
use  of  some  of  our  new  results.  These  examples  are  extensions  of  the  ones  given  pre- 
viously in  (37)  by  having  the  opening  range  of  battle  be  less  than  the  minimum  of  the 
two  maximum  effective  ranges  for  the  two  weapon  systems.  They  are  motivated  by  the 
work  of  S.  Bonder  [6],  (8),  [111  on  the  value  of  range  capabilities  and  mobility  for 
weapon  systems  in  combat  described  by  Lanchoster-typc  equations  of  modern  warfare. 
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The  modelling  context  of  these  examples  Is  that  of  weapon  systems  with  (a)  different 
range  dependencies  of  lethality  of  each  side's  fire  (but  the  same  maximum  effective 
range)  and  (b)  linear  attrition-rate  coefficients  but  different  effective  ranges. 

As  in  [37]  we  consider  a constant-speed  attack  on  a static  defensive  position 
with  the  combat  dynamics  described  by 

dx/dt  - -o(r)y  - -o0(l-r/Ra)my , dy/dt  - -0(r)x  - -0o(l-r/Rg)nx,  (45) 

where  R and  Ra  denote  the  maximum  effective  ranges  of  the  Y and  X weapon 
systems,  respectively  (i.e.  a(r)  - 0 for  r > Rq)  . Range  is  related  to  time  by 
r(t)  ■ Rq  - vt,  where  Rq  denotes  the  opening  range  of  battle.  Several  range  depend- 
encies for  an  attrition-rate  coefficient  are  shown  In  Figure  2,  and  an  opening  range 
less  than  the  weapon  system's  maximum  effective  range  is  indicated.  The  parameters 
of  the  attrition-rate  coefficients  a(r(t))  and  0(r(t))  in  (45)  are  readily  related 
to  those  of  a(t)  and  b(t)  In  (4)  and  (5 >.  For  example,  we  have  for  (4)  that 
kfi  - o0(v/Ra)m,  - 0o(v/Rg)n,  and  C - (Ra-R0)/v,  where  Rq  - Rp. 

Numerical  results  are  shown  In  Figures  3,  4,  and  5.  The  force-level  trajectories 
have  been  generated  by  a digital  computer  program  usln&  the  auxiliary  Lanchester  func- 
tions (i.e.  Fy  and  for  (31)  and  h,  w,  H,  and  W for  (44)).  These  functions 
are  particularly  convenient  for  such  digital  computer  work.  The  computer  routines 
were  checked  against  the  numerical  results  given  previously  in  [37].  In  this  work  we 

have  taken  the  opening  range  of  battle  Rn  < minimum  (R  ,R„).  Nunerlcal  values  for 

U OP 

battle  parameters  (except  those  for  Rn,  R , and  R.)  are  the  same  as  those  used 

U Cl  p 

for  the  examples  considered  in  [37]. 

For  Figures  3 and  4 both  weapon  systems  have  the  same  maximum  effective  range 
(i.e.  Rfl  ■ Rg) . As  done  for  the  plots  in  [37],  we  have  held  <»q  ■ a(r-O)  and  0q 
constant  and  have  varied  the  exponents  m and  n which  control  the  range  dependencies 
of  a(r)  and  8(r).  With  the  exception  of  Rq,  all  the  battle  parameters  for  the 
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force  weapon  system  kill  rate  at  zero  force  separation  (range).  3.  The  opening  range 
of  battle  is  denoted  as  R-  ■ 1250  meters  and  (as  shown)  R-,  < R .] 


X force 
level 
x(t) 


Y force 

level 
y(t ) 


15 
10 

5 

2 50  500  750  1000  1250 

Force  separotion  Jb  ( meters  ) 

250  500  750  1000  1250 

30 
20 
10 


Figure  3.  Force-level  trajectories  of  X and  Y forces  for  various  combina- 
tions of  the  exponents  m and  n in  the  power  attrition-rate 
coefficients  for  Rq  ■ 1230  meters  , Rq  ■ Rg  ■ 2000  meters, 

- 0.06  X casualties/(rainutes  ' Y unit)  , Bq  - 0.6  Y casualties/ 
(minutes  • X unit)  , v = 5 raph,  Xq  ■ 10,  and  ■ 30.  The  exponent 
combinations  are  denoted  as  m:n  in  the  figure,  and  the  symbol  * 
denotes  the  end  of  a force-level  trajectory  due  to  the  annihilation 
of  the  enemy  force. 
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X force 


Y feres  20 


Force  separation  Sv  (meters) 


2 50 


500 


750 


I00( 


Figure  4.  Force-level  trajectories  of  X and  Y forces  for  various  com- 
binations of  the  exponents  m and  n in  the  power  attrition-rate 
coefficients  for  the  same  parameter  values  chosen  for  Fig.  3 
except  that  ■ Rg  ■ 1500  meters.  The  symbol  conventions  are 
also  the  same  as  in  Fig.  3. 


X force 
level 
x(t) 


250  500  750  1000  1250 


Figure  5.  Force-level  trajectories  of  X and  Y forces  for  various 

effective  ranges  R^  of  the  X force  weapons  with  linear 

attrition-rate  coefficients  for  R - 1500  meters  and  the 

cj 

same  values  of  the  other  parameters  listed  in  the  legend 
of  Fig.  3.  The  symbol  * has  the  same  meaning  as  in  this 
figure. 
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curves  shown  in  Figure  3 are  Che  same  as  those  used  for  the  corresponding  example 

previously  considered  in  (37]  (i.e.  R ■ R ■ 2000  meters).  Consequently,  corre- 

o p 

sponding  force-level  trajectories  are  similar  with  greater  "separation"  shown  here 

between  curves.  In  Figure  3 with  the  opening  range  of  battle  Rq  ■ 1250  meters,  the 

curves  corresponding  fo  the  constant-coefficient  case  (i.e.  m • n ■ 0)  are  exactly 

the  same  (for  the  same  time  intervals)  as  those  shown  in  (37)  with  Rq  ■ 2000  meters. 

Other  battle  trajectories  with  m,n  > 0 decay  faster  in  Figure  3 than  they  did  in 

(37]  because  the  "intensity"  of  combat  is  greater  (i.e.  as  a function  of  time  the 

attrition-rate  coefficients  are  larger  here  than  they  were  in  (37]  (see  Figure  2)). 

As  noted  in  ( ? 7 ] , knowledge  of  the  range-  (or  time-)  dependence  of  weapon 

system  kill  capability  is  essential  for  forecasting  the  battle's  outcome  from  the 

initial  trend  of  battle.  For  example,  compare  the  outcomes  for  curves  denoted  as 

1:0,  1:1,  and  1:2.  G.  Clark  (13]  has  developed  nethodology  for  estimating  such 

capability  from  the  output  of  a (high  resolution)  Monte  Carlo  combat  simulation.  The 

"compounding"  nature  of  attrition  over  time  for  (45)  is  evident  from  the  curves  shown 

in  Figure  3:  a small  advantage  In  numbers  and  firepower  becomes  magnified  over  time. 

Similar  battle  curves  are  shown  in  Figure  4 for  the  same  parameter  values  except  that 

R - Rn  ■ 1500  meters.  Observing  that  for  m ^ 1 we  have  a(r;R  ) < a(r;R  ) » 

a p a a 

R < R , we  may  consider  that  the  combat  is  less  "intense"  for  such  pairs  of  battle 

a a 

trajectories  in  Figure  4 than  for  those  shown  In  Figure  3. 

In  Figure  5 wo  show  the  effect  of  increasing  the  effective  range  of  the  defender's 

weapons.  [The  X force  may  be  considered  to  be  the  defender.]  For  these  computations 

(using  (44)  and  the  auxiliary  Lanchesrer  functions)  we  have  held  the  opening  range  of 

battle  constant  at  Ry  = 1250  meters  and  the  maximum  effective  range  of  the  Y weapon 

constant  at  Rfl  B 1500  meters.  As  in  ( 37 ] , both  attrition-rate  coefficients  depend 

linearly  on  range  (i.e.  m = n B 1 in  (45)),  and  8^  have  been  held  constant, 

and  R„  has  been  varied.  As  shown  in  Figure  5,  we  quantitatively  see  the  benefit 
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Figure  6.  Classification  of  Lanchester-Type  Equations  of  "Modern  Varfare 
and  Their  Ease  of  Solution  by  Analytical  Methods 
(after  L.  Bertalanffy  [5]). 


dx/dt  « -a(t)y  - 0(t)x 


with  x(t*0)  - Xq, 


(49) 


dy/dt  - -b ( t ) x - a(t)y  with  y(t-O)  ■ yn, 


where  the  attrition-rate  coefficients  are,  of  course,  nonnegative.  In  this  case  the 
X force-level  equation  is  given  by 


d*x 


+ (o(t)+6(t)  - + <«<OKt)+»<t>  - ®(t)b(t))x  • 0.  (50) 


dr  6(0  ■ 


dt 


The  theory  of  this  paper  (also  that  of  (37])  may  be  applied  to  (49)  (or  equivalently 
(50)).  However,  to  us  this  appears  to  represent  essentially  the  limits  of  fiultful 
analytic  investigation.  Except  for  special  cases,  the  equations  for  combat  between 
heterogeneous  forces  do  not  appear  to  be  amenable  to  analytic  solution.  To  be  sure, 
one  can  write  down  "symbolic  solutions"  (such  as  the  matrix  exponential  for  constant 
attrition-rate  coefficients  (see  [11],  [37])),  but  these  are  apparently  not  computa- 
tionally useful. 


10.  Discussion. 

In  this  section  we  discuss  the  significance  and  applications  of  the  results  of 
this  paper.  These  results  may  be  used  to  facilitate  parametric  analysis  of  the  dynamic 
combat  interactions  between  two  homogeneous  forces  with  time-  (or  range-)  dependent 
weapon  systcr  capabilities.  Such  models  are  of  particular  interest  in  light  of  the 
work  of  S.  Bonder  [7],  (9]  and  others  [3],  [11],  [31],  [32]  on  the  prediction  of 
Lanchester  attrition-rate  coefficients  from  weapon  system  performance  data  and  the 
work  of  G.  Clark  [13]  on  the  estimation  of  such  (time-dependent)  coefficients  from 
Monte  Carlo  simulation  output.  Additionally,  our  new  theoretical  results  may  be  used 
to  solve  Grubbs  and  Shuford's  [19]  new  probabilistic  formulation  for  Lanchester  conbat 
theory  (see  equation  (3)).  A further  discussion  of  applications  is  to  be  found  in  [37] 
(see  also  [11]). 

We  have  presented  a general  mathematical  theory  for  the  solution  of  variable 
coefficient  Lanchester- type  equations  of  "modern  warfare"  for  combat  between  two 
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homogeneous  forces.  These  results  allow  one  to  extend  results  given  in  [37]  that 
applied  under  rather  restrictive  conditions  (e.g.  opening  range  of  battle  equal  to  the 
minimum  of  the  two  maximum  effective  ranges  of  the  weapon  systems)  (see  also  [11]). 

It  has  been  shown  that  in  general  the  deterministic  time  histories  of  the  X and  Y 
force  levels  (i.e.  x(t)  and  y(t))  may  be  expressed  in  terms  of  four  complementary 
Lanchester  functions  (and  not  two  as  in  the  constant-coefficient  .case) . The  mathe- 
matical properties  of  these  general  Lanchester  functions  that  facilitate  analytic 
solution  ha/e  been  given.  We  have  given  particular  attention  to  the  determination  of 
the  conditions  on  attrition-rate  coefficients  that  allow  relatively  simple  analytic 
results. 

Motivated  by  the  analytic  solutions  of  such  equations  for  certain  functional 
forms  of  attrition-rate  coefficients  of  interest  (i.e.  power  attrition-rate  coefficients 
and  linear  attrition-  rute  coefficients  with  different  effective  ranges),  we  have  pro- 
posed new  mathematical  functions  (and  given  their  mathematical  properties)  which  could 
be  tabulated  in  the  future.  This  would  allow  analysts  to  generate  numerical  results 
for  such  variable  attrition-rate  coefficient  combat  formulations  with  somewhat  the 
same  facility  as  one  can  for  the  constant-coef licient  case  and  thus  aid  in  parametric 
analyses.  This  means  that  numerical  results  for  a single  battle  can  be  used  to  generate 
numerical  solutions  for  an  entire  family  of  battles.  Such  results  for  the  reference 
battle  need  not  even  be  generated  by  analytic  means  but  could  be  developed,  for  example, 
by  numerical  integration.  Thus  one  can  combine  the  general  th  -oretical  results  given 
here  with  tabulations  generated  by  finite  difference  methods.  We  have  discussed  what 
appear  to  be  the  limitations  of  such  analytic  investigations  as  far  as  complexity  of 
combat  dynamics. 

The  results  of  the  paper  at  hand  and  [37]  hopefully  provide  the  theoretical 
foundations  for  the  generation  of  force-level  trajectories  via  "analytic"  solutions 
for  variable-coefficient  combat  formulations  (Lanchester-type  equations  of  "modern 
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warfare")  with  approximately  the  same  degree  of  thoroughness  in  these  foundations  as 
previously  existed  for  the  constant-coefficient  case.  In  other  words,  we  have  stressed 
computational  procedures  here.  For  qualitative  insights  into  the  dynamics  of  combat 
via  the  force  ratio  reference  should  be  made  to  the  recent  paper  by  Taylor  and  Parry 

(40) . 

11.  Notes. 

1.  This  is  because  of  the  complexity  of  obtaining  analytic  results  from  stochastic 
formulations.  The  computational  limitations  of  Monte  Carlo  simulation  (especially  high 
resolution  or  for  combat  between  large  units)  are  well  known.  The  work  of  G.  Clark 
shows  (see  pp.  102-103  of  [13])  the  complexity  of  an  analytic  solution  for  even  the 
most  idealized  combat  situation  modelled  as  a continuous  parameter  Markov  chain  (see 
[13]  for  further  references).  Moreover,  Bonder  and  Farrell  [11]  report  excellent 
agreement  between  simulation  results  and  those  for  a corresponding  deterministic 
Lanchester-type  model  (numerical  solution  generated  by  finite  difference  methods  (see 

[41] )).  In  this  sense  modem  high  speed  digital  computers  have  facilitated  the  applica- 
tion of  general  systems  theory  (see  [5],  especially  pp.  17-20)  to  military  systems. 

Of  course,  verification  of  such  models  (as  with  any  combat  model)  is  an  unresolved 
question  (see  [10];  further  references  are  given  in  [37]). 

2.  Although  finite-difference  mehtods  (see  [41])  and  a modem  high  speed  digital  com- 
puter can  generate  approximate  numerical  solutions  to  Lanchester-type  equations  with 
theoretically  any  degree  of  desired  accuracy,  it  nevertheless  is  of  interest  to  have 
analytic  solutions  available  (see  [28]),  if  for  no  other  reason  than  to  be  able  to  check 
the  adequacy  of  finite-difference  approximation.  Considering  the  ease  of  generating 
such  numerical  solutions,  it  is  curious  that  Barfoot  states  that  (p.  888  of  [3]) 
"existing  Lanchestcr  models  of  combat  have  assumed  constant  attrition-rate  coefficients 
for  the  sake  of  simplicity."  This  statement,  in  fact,  apparently  goes  back  (via 
Bonder  (p.  231  of  [7]))  to  Dolansky  (p.  345  of  [16]),  who  also  stressed  t lie  (at  that 
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time,  1964)  lack  of  valid  means  for  determining  such  coefficients.  As  one  of  the  authors 
has  shown  [34],  for  a constant  ratio  of  attrition-rate  coefficients  the  solution  to 
such  variable-coefficient  equations  for  combat  between  two  homogeneous  forces  is  no 
more  complicated  than  that  for  constant  coefficients. 

3.  The  influential  military  philosopher  of  the  19th  century,  Carl  von  Clausewitz 
(1780-1831)  stated  in  his  classic  work  On  War  (Vom  Krlege)  (p.  276  of  [14J),  "The 
best  Strategy  is  always  to  be  very  strong,  first  generally  then  at  the  decisive  point. 
...There  is  no  more  imperative  and  no  simpler  law  for  Strategy  than  to  keep  the  forces 
concentrated." 

4.  There  is,  moreover,  far  from  universal  agreement  as  to  what  are  the  significant 
variables  which  describe  the  combat  process  and  can  be  used  to  predict  its  outcome. 

For  some  other  views  see  [21]  and  [26]. 

3.  The  differential  equation  under  consideration  could  not  be  found  among  the  443 
linear  second  order  equations  tabulated  in  [24]  or  the  596  tabulated  in  [ 30  j . 

6.  For  a discussion  of  algebraic  addition  theorems  see  Chapter  II  of  [20].  Harris 

Hancock  gives  the  following  theorem  '(see  p.  37  of  [20]):  a necessary  and  sufficient 

condition  for  a single-valued  analytic  function  f(z)  to  possess  an  algebraic  addition 
theorem  is  that  there  exist  between  the  function  f(z)  and  its  first  derivative  an 
algebraic  equation  whose  coefficients  are  independent  of  the  argument  z.  No  more 
comprehensive  result  could  be  found  in  the  more  recent  monograph  by  Aczel  [2]. 

7.  A function  similar  to  Fy(x)  was  introduced  by  Ludwig  SchlMfli  in  1867  [33]  and 
another  appears  in  a posthumous  fragment  by  William  Kingdon  Clifford  (1845-1879)  (see 
pp.  346-348  of  [15]).  Greenhill  [18]  suggested  that  such  a function  would  be  conven- 
ient for  certain  engineering-like  problems,  although  he  was  severely  criticized  by 


35 


Watson  (see  p.  91  of  (42]).  Since  tabulations  of  none  of  the  above  or  modified  Bessel 
functions  of  the  first  kind  of  fractional  order  (except  for  the  restrictive  set  of 
values  v ■ ±1/4, ±1/3, ±1/2, ±2/3, ±3/4)  exist,  it  would  seem  appropriate  to  introduce 
the  LCS  functions  as  we  have  done. 
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